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Abstract 

Biased random walk has been studied extensively over the past decade especially in the transport and communication 
networks communities. The mean first passage time (MFPT) of a biased random walk is an important performance indicator 
in those domains. While the fundamental matrix approach gives precise solution to MFPT, the computation is expensive and 
the solution lacks interpretability. Other approaches based on the Mean Field Theory relate MFPT to the node degree alone. 
However, nodes with the same degree may have very different local weight distribution, which may result in vastly different 
MFPT. We derive an approximate bound to the MFPT of biased random walk with short relaxation time on complex network 
where the biases are controlled by arbitrarily assigned node weights. We show that the MFPT of a node in this general case 
is closely related to not only its node degree, but also its local weight distribution. The MFPTs obtained from computer 
simulations also agree with the new theoretical analysis. Our result enables fast estimation of MFPT, which is useful 
especially to differentiate between nodes that have very different local node weight distribution even though they share 
the same node degrees. 
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Introduction 

Scale-free node degree distribution, small network diameter, 
large clustering coefficients - these are common properties found 
to be present in complex networks arising from seemingly 
disparate fields such as biology, computer science, cosmology, 
etc. [1,2]. It is widely believed that there should be common 
underlying principles behind the formation of these networks that 
resulted in the observed properties. As such, complex networks 
have received much research attention during the past decade. 

One of the studies pertaining to complex networks is the 
network efficiency and capacity analysis [3-7]. In these studies, 
answers to the questions such as 'how fast can a message be 
delivered to a given destination' and 'how many packets may be 
generated in the system before congestion arises' are essential in 
understanding the performance of a network [5,7,8]. While there 
are existing results on routing and flow balancing in networks with 
certain topologies [9], they usually assumed the underlying 
topology or the knowledge of link formation mechanism. There 
are cases where the assumptions do not hold, for example, the 
animal foraging strategy [10] and the web searching process as 
depicted by the PageRank algorithm [11]. Approaches based on 
random walk can be applied when the detailed information of 
network formation mechanism is absent. Quantities such as 
stationary distribution and mean first passage time (MFPT) are important 
as they can be used to answer the questions about the performance 
of networks as mentioned above. 

The concept of random walk has also been applied in social 
networks. Even-Dar et al. [12] studied the process of spreading 
influences in social networks by means of Voter Model and 



showed that the pathways in which the influences propagate are 
equivalent to series of random walks. Thus, the MFPT to a node 
A, yields the expected time for the other nodes in the network to 
be influenced by node A. Selecting a node with low MFPT for 
spreading the news could result in fast information propagation. 

In fact, MFPT can be calculated by using the fundamental 
matrix method [13]. However, the computation involves multiple 
matrix multiplications. When the method is applied on large scale 
networks with millions of nodes, the computation becomes 
practically infeasible. Moreover, the solutions obtained from the 
fundamental matrix approach are too generic and hard to 
interpret. For instance, it is unclear which factors, be it node 
degree, eigenvalue, local connectivity, or others, govern the MFPT 
by just looking at the solution expression. Further research is 
needed to better characterize MFPT and to reduce the compu- 
tational cost. 

In [10], Condamin et al. showed a mean first passage time 
analysis using the pseudo Green function. They related MFPT to 
the network size and diameter. The general applicability of their 
result to non-fractal networks has been discussed in [14]. Fronczak 
et. al [8] applied the mean field theory to study the MFPT based 
on the Erdos-Renyi (ER) random graphs and networks generated 
by using the Barabasi-Albert (BA) preferential attachment model. 
Lau et. al [14] showed asymptotic analysis of the first passage time 
of unbiased random walk for a class of networks with short relaxation 
time by using the mean field theory. However, the solution relates 
MFPT to the node degree alone. There are cases where the nodes 
share the same node degree while having vastly different local 
topology such as those depicted in Figure 1. 
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Figure 1 . Nodes with same degree may have very different local connectivity. The figure shows two examples where a node with degree 3 

may be part of a sparsely-connected star network or a densely-connected clique. 

doi:10.1371/journal.pone.0093348.g001 



Random walks with short relaxation time are also known as 
random walks with non-compact exploration [15-17]. In [15], 
Benichou et al. presented the conditions for which a random 
walk on fractal network falls into the compact or non-compact 
exploration regime. In [16], Hwang et al. presented a fundamen- 
tal result on the relationship between the node degree and the 
MFPT of heterogeneous networks, i.e., networks with power-law 
degree distribution. They focused on uniform random walks 
without self-loops. It was shown that the general trend of MFPT 
with respect to the node degree exhibit cross-over behaviour that is 
governed by the spectral dimension d s and the exponent of the degree 
distribution y. For the case d s < 2, which is known as the compact 
exploration regime (see also [15]), the MFPT was found to be 
independent of node degree. For d s >2, the MFPT is related to the 
node degree via a power-law function where the exponent is 
determined by d s and y. 

While the focus of [16] is to show the general trend of MFPT 
with regard to node degree and spectral dimension, we focus on 
explaining the differences in MFPT for nodes with same degree 
when the underlying random walk falls into the non-compact 
exploration regime. We generalize the FPT analysis as discussed in 
[14] for a class of random walks with short relaxation time where 
the nodes have arbitrary weights. In addition, we analyse the first 
passage time at an improved level of precision by incorporating 
exact solutions to the stationary distribution. This enables a more 
detailed view of the neighbourhood around the target node. The 
new expression allows us to differentiate between the case where 
the neighbourhood of a node is sparsely connect and the case 
where it is densely connected as depicted in Figure 1 . We will show 
that the MFPT of a node is closely related to the local weight 
distribution around it. By changing the node weights, the spectral 
dimension and hence the relaxation time will be affected, we have 
also conducted simulations to study how the spectral dimension is 
affected by varying the weight assignments. They will be discussed 
in more detail in Discussions section. 

Analysis 

Outline 

In this section, we will show the detailed analysis for 
approximating the FPT decay rate and MFPT of random walks 
with short relaxation time. First, we will show the stationary 
distribution for random walk with arbitrary node weight 
assignment scheme. Then, we apply a flow-based heuristics to 
estimate the quasi-stationary distribution of the random walk when 



a sink node is introduced. Next, we show that the FPT distribution 
follows an exponential decay, and further show that the decay 
exponent is related to the quasi-stationary distribution of the 
neighbours of the target node. Finally, we obtain the MFPT by 
approximating the integration over the FPT distribution. 

Stationary and Quasi-stationary Distribution 

A complex network is modelled by a connected and undirected 
network G = (V,E,W), where V denotes the set of nodes with 
\V\=N, E denotes the set of edges (ij), and W denotes the 
weights assigned to the nodes. An edge (ij) represents the 
existence of relationship between nodes i and j. For the sake of 
simplifying expression, we assume without loss of generality that 
self-loops are present for all nodes, i.e., (ij)eE for /e[l,7V]. We 
define a biased random walk on G with arbitrary positive weights 
assigned on nodes. A transition from node i to a neighbouring 
node j is based on the following transition rule: 



Pr{i^j)- 



(1) 



where w z - denotes the weight of node /, and n(i) — {h\(i,h)eE} 
denotes the neighbourhood of node i. Pr(i—>j) denotes the 
probability of a random walker at node i moving to node j at the 
next time step. When self-loop is present, the staying probability of 
a node is not a constant but is dependent on the local weight 
distribution. Let the transition matrix be denoted by 
P=[Pr(i^>j)], and Pf(t) denote the probability of the random 
walker appearing at node i exactly at timestep t. The master 
equation is given by: 



jen(i) 



(2) 



Recall that a Markov Chain (MC) is said to be regular if 3q gZ + , 
V(iJ)p( q \i^>j)>Q where P iq \i^>j) denotes the (ij) entry of the q- 
th order transition matrix P q . The probability distribution of a 
regular MC will converge to a unique stationary distribution n 
regardless of the starting position as >oo. The relaxation time t of a 
Markov Chain is the time for the state probability distribution to 
be close to the stationary distribution, i.e., the standard deviation 
of P T is bounded by 1 je where e denotes the Euler Number (see 
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Figure 2. Plots of empirical first passage time distribution against theoretical prediction according to the approximate bound 
given by lneq.(13) for different networks and weighting factors. Each row corresponds to a network in the following order: Actor, BA, ER, 
and arXiv. The columns, from left to right, correspond to a = — 1, 0, 1 respectively. For most cases, the tail of the first passage time distribution can be 
predicted fairly accurately except for Figure 2(j), which is due to the high relaxation time as shown in Table 2. 
doi:1 0.1 371 /journal. pone.0093348.g002 



[18] Chapter 12). The relaxation time of a Markov Chain is given 
by i J; ? where hi is the second largest eigenvalue [18]. 

In the networks that we consider, since every node has a self- 
loop and the graph is strongly connected, the corresponding 
Markov Chain is regular. In the following, we show the exact form 



of stationary distribution for the biased random walk as defined in 
Eq.(l). 

Let Wi be the weight of node i, wj be the neighbourhood weight 
of node i given by w\ = ^iYj^n{i) w h and = ^J2M ■ The 
stationary distribution of the biased random walk is given by 
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Table 1. Summary of networks studied. 





Network 


# nodes 


# edges 


CC 


Dm 


Source 


arXiv 


4158 


15501 


0.5569 


17 


[21] 


BA 


4158 


31136 


0.0174 


5 


[1] 


ER 


4158 


21020 


0.0022 


7 


[24] 


Actor 


968 


13324 


0.6751 


9 


[1] 



C denotes the average clustering coefficients and Dm denotes the network diameter. 
doi:1 0.1 371/journal.pone.0093348.t001 



71/ = limP'(0 = 



Ex ™x E yen (x) W y 



N- (w f ) 



(3) 



The expression can be verified by applying the equilibrium 
condition on the master equation as given by Eq. (2). 

Let Vd denote the sink node and let G Vd denote the resultant 
graph from designating node Vd as sink node. Let <$>(t) denote the 
corresponding random walk probability distribution and <j) v (t) 
represent the probability for a random walker to be present at 
node v at time t. The following set of rules describes the new 
random walk on vJ v , '. 



where (/> to tal0) = £veK^v(0 denotes the total survival probability 
at time step t. 

We approximate the quasi-stationary distribution of G Vd around 
the sink node with a flow-based heuristic [14] described in the 
following. We begin with the stationary distribution n in G. Under 
the stationary distribution, based on Eq.(3), the probability of the 



random walker traversing an edge (w,v) is P n (u,v) = 



N(w') 



Next, 



we treat an undirected edge as a combination of in-link and out- 
link. We remove the out-links from Vd and hence making Vd a sink 
node. Thus, the equilibrium will be broken and nodes uen(vd) will 
have their 'flow' constantly drawn by Vd- Finally, we approximate 
the quasi-stationary distribution of such nodes by discounting the 
probability of utilizing the edge (u,Vd): 



0v(O) = <5v, 

^(0 = 0 



uen(v) 1^ W r 

ren(v) 



(initial condition) 
(sink node) 

(transition rule) 



(4) 



where Ve denotes the source node and d v „ = 1 for 



and 0 



otherwise. For connected network where Vd is reachable by every 
node, (j> v (t) tends to zero when t tends to infinity. However, for 
random walks with short relaxation time [14,19], the conditional 
probability distribution Pj(t\T>t) (i.e. conditioned on the survival 
of the random walker) will converge as t»T, where T denotes the 
time to reach the sink node. As such, the converged conditional 
probability distribution is called the quasi-stationary distribution 
denoted as *P V . For t»T, we have the following approximation: 



4 total® 



(5) 



i(l )n u = (l -Pr(u^v d ))n u 



(6) 



Asymptotic First Passage Time Analysis 

To obtain FPT, we re-designate the destination node Vd as a 
sink node such that the random walk process terminates once the 
random walker moves into the sink node. The time taken for it to 
be absorbed into the sink node is then the same as the FPT. Let 
F(vd,t\v s ) denote the probability of visiting node Vd at timestep t 
for the first time with v s as the starting node. Then by definition, 
we have 



F(VdA Vj) = hotalit - 1) - ^ total® 



(7) 



i.e. the first passage probability at time t is given by the difference 
in the total survival probability between time t and t — 1 . 

On the other hand, we can also obtain the first passage 
probability by using the transition rule. 



Table 2. Pearson's correlation coefficient for MFPT. 





Network 


Correlation 






Relaxation time r 










0 


1 


a=-\ 


0 


1 


arXiv 


0.2619 


0.7435 


0.9269 


555.56 


555.56 


1666.7 


BA 


0.9241 


0.9812 


0.9940 


3.5311 


2.2060 


2.0178 


ER 


0.4860 


0.9846 


0.9973 


16.502 


3.7327 


2.8417 


Actor 


0.0064 


0.9231 


0.9851 


476.19 


61.728 


26.041 



Overall, the correlations are high whenever the relaxation time is low. For the BA network, the relaxation time is consistently low and thus the correlation is extremely 
good. 

doi:1 0.1 371/journal.pone.0093348.t002 
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Figure 3. Comparison of MFPTs obtained from simulation against that obtained from using the bound shown in Ineq. (16). The 

column on the left compares the empirical MFPT to the results obtained by (i) using the bound in Ineq. (1 6); (ii) using the result presented in [8]. The 
column on the right shows the scatter plot and correlation between the empirical results and our proposed theoretical results. The rows correspond 
to the Actor, arXiv, and BA network respectively. While the result presented by Fronczak et al. [8] gives the general trend of MFPTs with respect to 
node degree, we find that the MFPTs for a given node degree are distributed across a wide range and cannot be fitted with a function of the node 
degree alone. The scatter plots show a strong correspondence between the empirical results and our proposed theoretical results. This is further 
supported by the high Pearson correlation coefficients which are shown on top of the scatter plots. 
doi:1 0.1 371 /journal. pone.0093348.g003 
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Table 3. Spectral dimension d s w.r.t. weighting factor a. 



d s 






arXiv 


BA 


ER 


-50 


0.042 


0.016 


0.026 


0.110 


-25 


0.068 


0.028 


0.142 


0.173 


-5 


0.360 


0.300 


3.760 


0.990 


-4 


0.476 


0.360 


4.500 


1.537 


-3 


0.590 


1.735 


5.066 


2.220 


-2 


2.000 


2.561 


00 


3.526 


— 1 


1 .932 


3.440 






0 


1.500 


2.505 


00 


00 


1 


1.456 


1.060 


CO 


CO 


2 


1.400 


0.667 


00 


00 


3 


1.025 


0.420 


00 


CO 


4 


0.920 


0.229 


00 


00 


5 


0.880 


0.220 


00 


4.378 


25 


0.180 


0.295 


00 


0.536 


-50 


0.107 


0.090 


00 


0.220 



The data is obtained by conducting simulation and fitting the exponent of RTO probability according to the definition given in Eq.(17). We mainly consider the range [— 
5,5] as they are mostly considered in the literature. The entries {-50, -25, 25, 50} are used to examine the effect on d s for large a. The entries with d s = oo are obtained 
by the observation that the RTO probability stabilized fairly quickly in less than 10 time steps. The data in this table is plotted in Figure 4. 
doi:1 0.1 371/journal.pone.0093348.t003 




Figure 4. Relationship between spectral dimension d s and the weighting factor a. The data used for drawing this figure is tabulated in 
Table 3. The lines corresponding to BA and ER network appear disconnected as they have d s = co for certain values of a and cannot be adequately 
shown in the figure. The spectral dimension generally peaks in the interval [-1,1] and drops significantly for a of greater magnitude. This is especially 
true for the BA and ER networks (from infinity to a finite value). In the extreme case, by setting a = oo, the 'random walk' is no longer random as the 
node with largest degree will always be chosen at every step. Similar reasoning also applies for oc = — go. Therefore, towards both extremes, we would 
expect the random walk to become more localized and hence falls into the compact exploration regime. 
doi:10.1371/journal.pone.0093348.g004 
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Figure 5. MFPT of the Actor network when a = 0. Even though the spectral dimension d s is less than 2, the MFPT is not found to be independent 
of node degree. Instead, the MFPT exhibits a power-law relationship with respect to the node degree. The disparity arises probably as a result of self- 
loops, which affects both the RTO probability distribution and the estimated d s . 
doi:10.1371/journal.pone.0093348.g005 



F(v d ,t\v s )= Ut-V ^ WVd 

uen(v d ) 2^ren(u) W r 



(8) 



Substituting the result back to Approx.(9) yields: 



F(v d ,t\v s ) a h d e- h d l '- 1) 



(13) 



For t<T, the first passage probability is dependent on the 
source-sink distance. Here, we focus on the asymptotic behaviour 
of the random walker for t»T. As such, we can apply Approx.(5) 
and after simplifying, we get: 



F(v d ,t\v s )^ lotal (t-l) Y, Vu ™ Vd 

uen(v d ) 2^ren(u) W r 



(9) 



Let p Vd - J2ua(v d ) " ■ 

Z— / ren(u) 



By combining Eq.(7) and 



Approx.(9), after simplification, we obtain the following recursive 
approximation: 



<l>total(i)~(l-Pv d )<l>total(t-l) 



(10) 



By definition, (j> tota i(0) = 1 as the random walker has not started 
moving and thus the survival probability is 1 . Solving the recursive 
formula yields: 



Thus, the first passage probability is approximately bounded 
below by an exponential function for t»T and the decay rate is 
given by fi vr To calculate the decay rate, we apply Approx.(6): 



v d 



uen{y d ) 



ren(u) W r J 



I UGll(v d ) \ 

E w "~ w 'd E y- - w 

I \uen(v d ) uen(v d ) ^ren(u) 



N(w') 



Wv d Euen(v d ) Pr(u^u) 



E, 



uen(v d ) 



-K Vd \ \-Pr(v d ^v d ) ^ Pr (u^u) 

\ uen{v d ) j 



(14) 



U')*(i-y 



<exp(-^0 



(11) 



(12) 



where the last step is just a simplification of the expression by using 
the definition as shown in Eq.(l) as we assumed earlier that every 
node in the network has a self-loop. Nevertheless, for the cases 
where self-loop is absent, the expression Pr(v— »v) can be 
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substituted for 



as self-loop played no role in our 



y ,,v 

derivation except for the regularity of the random walk. 

The result in Approx.(14) suggests that the decay rate of the first 
passage time distribution depends mainly on the sink node's 
stationary distribution and the transition probabilities around the 
sink node. 

Mean First Passage Time 

After obtaining the first passage time distribution, we can 
estimate the mean first passage time by using the approximation: 



T) = f2'F(v d ,t\v s ) 

1 

) 

tF(vd,t\v s )dt 

POO 

tF(v d ,t\v s )dt+ tF(v d ,t\v s )dt 



(15) 



for some c»t. In Approx.(15), the integral is separated into two 
parts - from 1 to c and from c to go for short-term and long-term 
behaviour respectively. While the first passage time analysis is 
based on the assumption t»T, our result is an upper bound on the 
first passage probability for small t. This is justified by the 
observation that the first passage time probability starts from zero 
initially, increases to a peak, then decreases exponentially (see 
Figure 2 and [15] for non-compact case) with the decay rate as 
presented previously. Thus, Ineq.(13) is an overestimation of the 
first passage probability for small t. Hence, we provide a lower 
bound of mean first passage time as follows: 



<r>> 



(1-//L>- / V'- I) A+ e-W^dt 
l d Ji 



(16) 



= 1 + 



v d 



Although Ineq.(16) is a lower bound of the mean first passage time, 
results from simulations show that it is strongly correlated to the 
actual MFPT. 

In order to calculate the MFPT, an edge will be followed exactly 
twice to obtain the weight of the connected neighbour. Thus, 
assuming the network is represented in adjacency list format, the 
time complexity of the approximation is given by 0(N + M). The 
time complexity to convert a network from adjacency matrix 
representation into adjacency list representation is bounded by 
0(N 2 ). Hence, the overall time complexity of the approach is 
bounded by 0(N 2 +M) for adjacency matrix representation. 
Nevertheless, in terms of computational efficiency, the approach is 
a great improvement from the fundamental matrix method [13] 
which has time complexity of 0(N 3 ) or 0(7V 2 " 376 ) depending on 
the actual implementation of matrix operations [20] . Thus, our 
approach can be used to estimate MFPTs of the network quickly 
especially for large scale networks. 

Results 

We verify the theoretical results on two real world networks and 
two artificially generated networks: (i)arXiv General Relativity and 



Quantum Cosmology collaboration network, obtained from the 
Stanford SNAP website [21], (ii) Barabasi- Albert (BA) preferential 
attachment network [1,22], (iii) Erdos-Renyi (ER) random graph, 
and (iv) Actor collaboration network from the Barabasi lab [1]. 
Because the analysis only applies to connected graphs, we used the 
largest component of the arXiv network. The two generated 
networks are chosen to test the theoretical results on networks with 
different structure and edge density. The ER network is generated 
with p = 0.0011, while the BA network is generated with 
m = No = l. The Actor network is constructed from the first 
1 000 records of the database, where each record consists of the 
actors who collaborated in the same movie. Table 1 summarizes 
the networks that we have examined. 

We apply the node weight assignment scheme given by Wi = kf, 
where k\ is degree of node i and a is an integer. This node weight 
assignment scheme is mainly studied in the network traffic 
community such as [5,7]. For our experiments, we mostly consider 
the range of a in [—1,1] except for the Actor network as other 
values with greater magnitude will result in random walks with 
exceedingly long walk lengths. 

The experiments are conducted as described below. Firstly, for 
a given network, we choose 20 nodes randomly as the source 
nodes. The weighting factor a is then fixed and the node weights 
are computed accordingly. For each source node, 250 times of full 
random walk simulations are conducted independently. A full 
random walk simulation starts with the random walker at the 
source node and terminates when every other node has been 
visited at least once. The first passage time to each node is 
recorded. The procedure is then repeated for other values of a. 

To obtain the first passage time distribution, we applied 
Gaussian Kernel Density Estimation on the first passage time 
statistics collected from the simulations. We randomly selected 
different source-sink pairs for each network and different values of 
a, and the results are plotted in Figure 2. As shown in the figure, 
for most cases, we can predict the tail of the first passage time 
distribution fairly accurately except for the case when a = — 1 in 
arXiv network. With reference to Table 2, we find that for 
a = — 1 , the relaxation time for the corresponding random walk is 
very high, and thus for certain sink nodes with high absorption 
rates, our theoretical result may not be applicable. 

In Figure 3, we compare the empirical MFPT to that predicted 
by Ineq. (16). Since we relaxed several expressions during the 
derivation, the values predicted by Ineq. (16) may not be of the 
same scale as the empirical result. Therefore, we renormalized the 
predicted values (P-set) with respect to the empirical result (E-set). 
The renormalization scheme is described as follows. First, we sort 
both the P-set and E-set in non-increasing order. Then we rescale 
the middle 90% of the P-set with respect to that of the E-set, i.e. 
obtain the rescale parameters (shearing and scaling) by ignoring 
both the upper and lower 5 percentile of both sets. Finally, the 
renormalization is applied to the whole P-set. 

As shown in Figure 3, we can observe that while the result by 
Fronczak et al. predicted the general trend of MFPT with respect 
to node degree, our result further refined the predicted values by 
examining local weight distribution. By zooming into a greater 
level of detail, our approach has revealed a useful relationship 
between local connectivity and the MFPT especially as highlighted 
in Figure 3(a),(c), and (e). The weighting factor a controls whether 
high degree nodes should receive greater attention or vice-versa. 
Thus, for small a, high degree nodes should be reached less often 
and hence greater MFPT. Surprisingly, even for a = — 2, there are 
cases where high degree nodes can be reached fairly quickly as 
depicted in Figure 3(a). By examining the structure of the Actor 
network, we find that many high degree nodes are connected with 
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bridge nodes of low degree. On the other hand, the MFPT of the 
BA network does not show this kind of pattern. We believe that 
this is due to the fact that the BA network possesses a much simpler 
structure than the Actor network. 

To investigate the overall performance of the bound given in 
Ineq. (16), we also calculated Pearson's correlation coefficient, and 
the results are summarized in Table 2. Overall, the correlations 
are high whenever the relaxation time is low. For the BA network, 
the relaxation time is consistently low and thus the correlation is 
extremely good. 

Discussions 

In [15,16], it was found that the random walk exploration can 
be divided into compact and non-compact regimes based on their 
spectral dimension. For the compact case, the random walker 
spends longer time travelling around certain neighbourhood and 
hence the relaxation time is long, while the non-compact case is 
the other way round. In our study, we focus on the non-compact 
case, i.e., random walk with short relaxation time. While the 
structure plays a part in deciding the dynamics of random walk, 
the weighting factor a also affects the mixing rates. To examine 
how much the assumption of short relaxation time holds when the 
weighting factor a is changed, we also investigated the relationship 
between a and the spectral dimension of the corresponding 
random walk d s . We numerically estimate d s by using the 
following relation [17]: 

/»<,(/)- r*' 2 (17) 

where P 0 (t) denotes the Return-to-origin (RTO) probability at 
time t [17]. To obtain the RTO probabilities, we conducted 
simulations with 100,000 random walkers which are placed 
randomly at / = 0. Each random walker walks independently for 
100,000 time steps, and the fraction of random walker returning to 
their respective starting node is recorded at each time step. Table 3 
and Figure 4 summarizes the result. 

As shown in Figure 4, the spectral dimension generally peaks in 
the interval [—1, 1], which explains the applicability of our result. 
While for the case of uniform random walk, it was shown that d s of 
both BA networks and random graphs are infinity [23], we found 
finite d s for a 7^0. In the extreme case, by setting a = — oc, the 
'random walk' is no longer random as the node with smallest 
degree will always be chosen at every step. For this scenario, the 
network will be broken into cycles where leaf nodes form the 
smallest cycle. Similar phenomenon also applies for a=oo. 
Therefore, towards both extremes, we would expect the random 
walk to become more localized and hence falls into the compact 
exploration regime where our result may not be applicable. 
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